Gaps between jets in double-Pomeron-exchange processes at the LHC 
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The possibility to measure jet-gap-jet final states in double-Pomeron-exchange events at the LHC 
is presented. In the context of the ATLAS experiment with additional forward physics detectors, 
cross sections for different experimental settings and gap definitions are estimated. This is done in 
the framework of the Forward Physics Monte Carlo interfaced with a perturbative QCD model that 
successfully reproduces standard jet-gap-jet cross sections at the Tevatron. The extrapolation to 
LHC energies follows from the Balitsky-Fadin-Kuraev-Lipatov dynamics, implemented in the model 
at next-to-leading logarithmic accuracy. 



o 

CsJ '. I. INTRODUCTION 

D . In a hadron-hadron collision, a jet-gap-jet event features a large rapidity gap with a high-pr jet on each side {pt^ 
' ^QCd)- Across the gap, the object exchanged in the channel is color singlet and carries a large momentum transfer. 
\^ , When the rapidity gap is sufficiently large, which requires a large collision energy {y/s ^ pr), the perturbative QCD 
P\| ■ description of jet-gap-jet events boils down to the exchange of a Balitsky-Fadin-Kuraev-Lipatov (BFKL) Pomcron 
[U, [l]. The first measurements carried out at the Tevatron [1,0] allowed preliminary tests of that QCD description, 
which were quite successful [11,13: paving the way for more thorough investigations at the LHC. 

In |5||: the full next-to- leading logarithm (NLL) BFKL kernel was implemented including all conformal spins and 
I . necessary coUinear improvements. In [^, the resulting parton- level hard cross section was interfaced with the HERWIG 
Monte Carlo [J, as was first proposed in ^ in the case of the leading- logarithm (LL) BFKL calculation, in order to 
f-H , take into account parton showering, hadronization, and jet reconstruction effects. The resulting model was able to 
reproduce the CDF and DO data, and to make predictions for the LHC, where many such events are expected to be 
' produced. 

' Unfortunately, due to the high luminosity and large pile-up environment, the jet-gap-jet measurement is hard to be 
' reproduced at the LHC. Instead, the so-called jet-veto measurement was proposed '9\, it amounts to vetoing on the 
I presence of additional jets in the inter-jet rapidity range, instead of strictly imposing a rapidity gap. However, doing 
' so significantly reduces the sensitivity to the BFKL dynamics [13, [HI- In this letter, we propose to take advantage 
fSI of the future ATLAS Forward Proton (AFP) detector fl^, in order to realize a true jet-gap-jet measurement. In 
this context, one would look for jet-gap-jet final states in double-Pomeron-exchange (DPE) processes, with both 
[ protons intact. In such events, rapidity gaps are easier to identify. We define a gap to be a rapidity interval devoid 
p\| ' of tracks with transverse momentum bigger than 200 MeV, for this reason we shall restrict ourself to central gaps, 
T-H within the ATLAS tracker acceptance. The modeling of such events can be obtained by implementing the parton- level 
h-jl . NLL-BFKL cross section into the Forward Physics Monte Carlo (FPMC) program |13l |. 

Compared to the original measurement, using jet-gap-jet final states in DPE processes implies cleaner events, not 
polluted by proton remnants, since those remain intact. This, in turn, allows to study jets with a larger rapidity 
difference, for which BFKL effects are supposed to be more important. On the theory side, the fraction of jet-gap-jet 
to inclusive di-jets events in DPE processes is larger than the corresponding fraction in non-diffractive processes. 
Indeed, in the former case the gap survival probability factor applies to both the jet-gap-jet and inclusive di-jets 
events, while in the latter case, it applies only to the jet-gap-jet events. 

The plan of the letter is as follows. In section II, we recall the phenomcnological NLL-BFKL formulation of the 
parton-level hard cross section, and we explain how it is embedded into the FPMC program. In section HI, we detail 
the experimental environment for the measurement, provided by the ATLAS experiment and its future forward proton 
detectors. Our predictions for the DPE jet-gap-jet cross section at the LHC are presented in section IV. Section V is 
devoted to conclusions and outlook. 
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FIG. 1. Left: diagram of a DPE event, with a diffractive final state made of two high-pT jets surrounding a large rapidity gap. 
Arjj is the rapidity interval between the jets whereas Ar^g denotes the rapidity gap. Right: the scheme of a DPE jet-gap-jet 
event in pseudo-rapidity plane. 



II. DOUBLE-POMERON-EXCHANGE JET-GAP-JET EVENTS IN FPMC 



To simulate DPE events, the FPMC program is used. This event generator is designed to simulate central particle 
production with two intact protons, and a hard scale present in the event. In this paper, we deal with the particular 
case when the centrally produced diffractive final state is made of two high-p^ jets {pt^ ^qcd) surrounding a large 
rapidity gap, as pictured in Fig. [1] 

Practically, in order to implement this specific final state in the FPMC, the function responsible for the calculation 
of the elementary parton-parton scatterings was modified. The 2-to-2 leading-order matrix elements ordinarily used 
for the case of central di-jet production, are replaced by 2-to-2 parton-level processes in the BFKL framework. This 
modification accounts for the i— channel exchange of a color-singlet gluon ladder with large momentum transfer, and 
therefore for the presence of a rapidity gap in between the final-state jets. Of course the collision energy ^/s should 
be big (-y/s ^ pt) in order for jets to be produced along with a large rapidity gap. As detailed below, the leading 
logarithm (LL) or next-to-leading logarithm (NLL) BFKL kernels are implemented, including all conformal spins and 
necessary coUinear improvements. 



A. 2-to-2 parton-Ievel processes in the BFKL framework 

In the gluon-gluon channel, the parton-level 2-to-2 hard cross section is given by 



dpi, IGtt 



MA7j,p'^)r , (1) 



where A{Ari,p'^) is the gg — > gg scattering amplitude. The incoming gluons, of longitudinal momentum fractions Xi 
and X2 respectively, carry zero transverse momentum, and therefore the final state partons have opposite transverse 
momenta: pt and — Pt- In addition, the rapidity gap coincides with the rapidity interval Ar/ = hi{xiX2S / p^) between 
the outgoing partons that will initiate the jets. The hadronization of the partons into jets produces a momentum 
imbalance and reduces the size of the rapidity gap, as is discussed in the next subsection. 

Let us now consider the high-energy limit in which the rapidity gap Arj is assumed to be very large. The BFKL 
framework allows to compute the gg — > gg amplitude in this regime 0, 0] : 

2 16iVe7ra|(p|) ^ f d7 [p^ - (7 - 1/2)"] exp {a(p|,)xe//[2p, 7, ^bDlAT?} 
^ '^'^^^ Cfp't ^£^^7 [(7 -1/2)2 -(p- 1/2)2] [(7 -1/2)2 -(p+ 1/2)2] ^ 
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with the complex integral running along the imaginary axis from 1/2— loo to 1/2+icx), and with only even conformal 
spins contributing to the sum. The running coupling is given by 

.(,1) = . [Mog (,^/A|,,)] H^^^ . (3) 

It is important to stress that in order to obtain formula ([2]), leading-order (non-forward) impact factors were used. 
In the case that one incoming particle is a quark instead of a gluon, the corresponding amplitude is obtained by 
multiplying the right-hand side of formula ^ by Cp/Nc- Note also that the Mueller- Tang prescription [3, [l^ was 
used to couple the t-channel BFKL Pomeron to the incoming colored partons. More details on these parton-level 
computations, such as the importance of the different conformal spins, or the size of the uncertainties due to the 
choice of the renormalization scale, can be found in jsj . 

In ©, the effective kernel Xeff{p,1,Ci) takes into account NLL-BFKL effects. For p — 0, the scheme-dependent 
NLL-BFKL kernels provided by the regularisation procedure XNLLil,^) depend on cj, the Mcllin variable conjugate 
to exp(Ar/). In each case, the NLL kernels obey a consistency condition ,16] which allows to reformulate the problem 
in terms of Xe//(7,Q;). The effective kernel Ye f f (7. is obtained from the NLL kernel xnll{i,^) by solving the 
implicit equation Xeff = XNLL{l:Ct: Xeff)- In [l2lill8|j the regularisation procedure has been extended to non-zero 
conformal spins and the kernel xnll(p,^, i^) was obtained from the results of [l^. The formulae needed to compute 
it can be found in the appendix of 17] o Then the effective kernels Xeff{p,1,0i) are obtained from the NLL kernel 
by solving the implicit equation: 

Xeff = XnlUp, 7, a Xeff) ■ (4) 

We recall that, in the LL-BFKL case, the formula for the jet-gap-jet cross section is formally the same as the NLL 
one, with the following substitutions in 

Xe//(p,7,a) ^ XLi(p,7) = 2?/'(l) - V'^^l - 7+ - V'^7+ , a ^ a = const, parameter, (5) 

where ■!/'(7) = dlogV {'-^) / d'y is the logarithmic derivative of the Gamma function. In this case, the coupling a is a priori 
a parameter, but the preferred value in order to fit the forward-jet data from HERA is 0.16 [l^]. This unphysically 
small value of the coupling is indicative of the slower Bjorken-a; dependence of the forward-jet data compared to the 
LL-BFKL cross section, when used with a reasonable a value. And in fact, the value a = 0.16 mimics the slower 
energy dependence of NLL-BFKL cross section (in this case the average value of a is about 0.25), which in the 
forward-jet case is consistent with data. 

B. Implementation of the NLL-BFKL formula in FPMC 

The parton-level calculation presented in the previous subsection leads by definition to a gap size equal to the 
interval in rapidity between the partons that initiate the jets. At the particle level, it is no longer true. Due to 
QCD radiation and hadronization, the jets have a finite size, and the gap size A?]^ is smaller than the difference in 
rapidity between the two jets Aryj (see Fig. [Ij. The implementation of the NLL-BFKL cross section ([1]) into the 
FPMC Monte Carlo takes these effects into account. Practically, we modified the function which implements the 
matrix element squared for elementary parton-parton scatteringj. Formula ([T]), which gives the BFKL da/dp^ cross 
section is too complicated to be implemented directly in FPMC since it involves an integration in the complex plane 
over 7, and it would take too much computing time to generate a reasonable number of events. To avoid this issue, 
we parameterized da/dp^ as a function of the parton pT and Ar] between both partons at generator level. Denoting 
^(Pt) ~ Q^(?'t)^'7/2, the parametrization used is 

^ = "fjy^ [a + bpT + +{d + epT + fy^) x z + {g + hpr) x + {i + j^) x + exp(fc + Iz)] . (6) 

This formula is purely phenomenological, not motivated by theory, and was just introduced to obtain a very good 
X^ while fitting ^ to the full expression of da/dp^. The fit was performed with 12 free parameters a — I, and 2330 
points were used. Formula dU) was implemented into FPMC with the fitted parameters. 



^ In the present study we shall use the S4 scheme (scheme 4 of Ref. fw^) in which Xnll is supplemented by an explicit a dependence. 

The results in the case of the S3 scheme are similar. 
2 This function is named HWHSNM in FPMC, as it was in HERWIG 
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FIG. 2. Scheme of the measurement concept: jets are registered in the central detectors, and protons in the very forward ones. 



III. EXPERIMENTAL ENVIRONMENT 

The ATLAS detector is located at the LHC Interaction Point 1 (IPl). It has been designed as a general purpose 
detector with a large acceptance in pseudorapidity, full azimuthal angle coverage, a good charged particle momentum 
resolution and a good electromagnetic calorimetry completed by a full-coverage hadronic calorimetry [2lj . The ATLAS 
tracking detector provides a measurement of charged particles momenta in the I77I < 2.5 region and the calorimeter 
covers I77I < 4.9. 

A crucial element of the DPE jet-gap-jet measurement is the possibility to tag the forward protons with the 
AFP detectors, see Fig. ^ The AFP detectors are being designed for proton-proton central diffractive or exclusive 
production measurements [l^. These detectors will be placed at 204 m and 212 m away from the IPl inside the 
" Hamburg beam-pipes" . These are special devices that allow to put detectors close to the beam and to control the 
distance between their edge and the proton beam. 



A. Forward Protons 

Since there are several LHC magnets between the IPl and the AFP detectors the proton trajectory depends not 
only on the scattering angle but also on the proton energy. Obviously, not all forward protons can be measured in the 
AFP detectors. Such protons can be either too close to the beam to be detected or can hit one of the LHC elements 
(a collimator, the beam pipe) before it reaches the AFP detector. This is included in Fig. [3l where the geometric 
acceptance of the AFP detector is shown. In this calculation the following factors were taken into account: the beam 
properties at the IP, the beam pipe geometry and the distance between the detector edge and the beam center. As 
can be observed, the region of acceptance is approximately limited by 0.012 < ^ < 0.14 and < 4 GeV, where 
^ = (1 — -E/-Eboam) is the relative energy loss and pt is the proton transverse momentum. 

Since both protons need to be tagged in the AFP stations, not all events can be recorded. In fact, the visible cross 
section depends on the distance between the AFP active detector edge and the beam centre (it will be varied during 
the runs according to the beam conditions). This dependence is presented in Fig. [3] (right). For the rest of the analysis 
a distance of 3.5 mm is assumed, which results in a visible cross section of about 1 nb (for a leading jet px > 40 GeV). 



B. Central Diffractive Jets 

The jets produced in the process under consideration will be measured in the central detector. In the following 
analysis, the cone algorithm with R = 0.7 was used for their reconstruction. To fulfill the ATLAS detector trigger, 
the leading jet is requested to have a transverse momentum greater than 40 GeV. This rather low value is realistic 
for the low pile-up LHC runs required (due to the gap reconstruction) to make this measurement. For reconstruction 
issues, the transverse momentum of the second leading jet is required to be greater than 20 GeV. In this analysis the 
two leading jets are required to be in the opposite pseudo-rapidity hemispheres and the rapidity gap is required to be 
symmetric around zero, i.e. from —\r]g\ to \rig\ with 1%! = Arig/2. 
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FIG. 3. Left: geometrical acceptance of the AFP detector as a function of the proton relative energy loss, ^, and its transverse 
momentum (pr)- Right: visible cross-section as a function of the distance between the detector and the beam centre (for 
leading jet with pr > 40 GeV). A distance of 3.5 mm would be at 15 sigma from the beam, this is the (standard) value that 
we choose. 



The main background for the DPE jet-gap-jet production will be DPE inclusive jet production. Indeed, in such 
processes, a gap between the jets can appear from fluctuations, but this background is significantly reduced by 
requiring large enough gap sizes (c/. Fig [J] (left)). Obviously, the bigger the gap size, the larger the DPE jet-gap-jet 
contribution. However, the cross-section falls steeply. The jet transverse momentum distribution for three different 
gap sizes {\r]g\ — 0.6, 1.0, 1.4) is presented in Fig. H] (right). Assuming both protons to be tagged in AFP, the good 
balance between signal to background ratio and visible cross-section is for the gap size |77g| > 0.5. 
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FIG. 4. Left: the gap size distribution for DPE jets and DPE jet-gap-jet events with and without the AFP tag requirement. 
We notice that for large enough gaps Arjg > 1, the gap-between-jets events are not dominated by fluctuations in di-jets events. 
Right: The jet transverse momentum distribution for different gap sizes. 
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FIG. 5. Predictions for the DPE jet-gap-jet to DPE jet cross section ratio at the LHC, as a function of the leading jet transverse 
momentum pr (left), and of the rapidity difference between the two leading jets Ar^j (right). For both plots, we assumed an 
integrated luminosity of 300 pb~^. 



IV. TEST OF THE BFKL MODEL AT THE LHC 



The DPE jet-gap-jet event ratio is defined as the ratio of the cross section for the DPE jet-gap-jet (JGJ) production 
to the DPE inclusive di-jet (Jets) production 



^ ajDPE JGJ) 

<j{DPE Jets) ' ^ ^ 

and is plotted in Fig.[5]as a function of the transverse momenta of the first-leading jet and as a function of the rapidity 
difference between the two leading jets, Ar/j. To verify the power of this test, the statistical errors corresponding to 
300 pb^^ of integrated luminosity were plotted. 

To take into account NLO QCD effects, absent in the FPMC program, the LO ratio obtained was corrected by the 
cross section ratio a{DPE LO Jet + +)/a{DPE NLO Jet + +), obtained by the NLO Jet-|-|- program This 
program, providing perturbative parton-level cross-sections was adapted for the diffractive applications as described 
in 23i (for di-jet production in p-l-p collisions, this program was independently validated by comparison with Frixione 
NLO [Ijl)- a crosscheck of the implementation of the resolved Pomeron model into NLO Jet-|--|-, the diffractive 
di-jet cross sections have been successfully compared to the predictions of the Monte Carlo generator RAPGAP [25l |. 

For the calculations, diffractive parton distribution functions (DPDF) from the HI 2006 Fit B were used [l^ (this 
fit arose from e-|-p HERA data on inclusive diffraction). The renormalization and factorization scales were set to be 
equal and identified to the leading jet transverse energy, i.e /ifl = = the jets were reconstructed by the cone 

algorithm with R=0.7 [13) i-e- by the same method as in FPMC, the number of fiavor was fixed to 5, and the strong 
coupling constant is taken at the 2-loop level with as{Mz) = 0.118, corresponding to that used in the evolution of 
the HERA parton densities. Finally, in the NLO Jet ++ calculations, a rapidity gap survival probability of 0.03 was 
applied for the LHC. It cancels in the ratio a{DPE LO Jet + +)/a{DPE NLO Jet + +). 

As far as the gap fraction a{DPE JGJ)/a{DPE Jets) is concerned, note that we did not consider an additional 
suppression factor for DPE jet-gap-jet production, on top of the 0.03 of DPE inclusive jet production. Therefore, in 
the predictions of Fig. [SJ all the rapidity gap survival probability cancel. We would like to point out that this last 
point is an assumption, DPE jet-gap-jet production could be subject to a bigger suppression than DPE inclusive jet 
production, due to the extra color-singlet exchange in the hard cross section. However, we do not expect this potential 
additional factor (on top of 0.03) to be large, due to the fact that extra soft interactions with the BFKL Pomeron are 
unlikely: the 2-to-2 hard scattering takes place on much shorter time scale compared to the soft interactions filling 
the rapidity gaps. In any case, this will be checked at the LHC and if necessary, our numbers can then be adjusted 
accordingly. 
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V. CONCLUSIONS 

In order to simulate the production of two high-px jets around a large rapidity gap in DPE processes in hadron- 
hadron collisions, we have embedded the parton-level BFKL calculations of |.5j into the FPMC program. The resum- 
mations of leading and next-to-leading logarithms are taken into account, and implemented through a renormalization- 
group improved kernel in the S4 scheme. The Mueller- Tang prescription is used to couple the BFKL Pomeron to 
colored partons, described with only leading-order impact factors. The implementation of that same NLL-BFKL cal- 
culation into the HERWIG Monte Carlo program led to a good description of all Tevatron data on standard jct-gap-jet 
production which provides a good foundation for this model. 

Due to the large pile-up environment, the jet-gap-jet measurement was not reproduced at the LHC, and the model 
predictions could not be checked. However the measurement can be done in DPE processes, which provides cleaner 
events not polluted by proton remnants, and consequently also gives access to larger di-jets with a larger rapidity 
difference, for which BFKL effects are more important. In addition, the fraction of jet-gap-jet to inclusive di-jets 
events in DPE processes is larger than the corresponding fraction in non-diffractive processes, since in DPE events 
one is not penalized by the gap survival probability, which applies to both the jet-gap-jet and inclusive di-jet cross 
section. 

After incorporating the parton-level BFKL cross sections into the FPMC generator, we obtained hadron-level 
results for jet-gap-jet final states in DPE events. In the context of the ATLAS experiment with additional forward 
physics detectors, we presented our cross section predictions for the LHC, for different experimental settings and gap 
definitions. They show that there will be enough statistics in 300 pb~^ of data to collect DPE jet-gap-jet events. This 
would provide an additional test of the BFKL Pomeron. 
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